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Abstract 


We determined the structures of silicon clusters in the 11-14 atom size 
range using the tight-binding molecular dynamics method. These calculations 
reveal that Sin is an icosahedron with one missing cap, Sii 2 is a complete 
icosahedron, Sii 3 is a surface capped icosahedron, and Sin is a 4-4-4 layer 
structure with two caps. The characteristic feature of these clusters is that 
they are all surface. 
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The size dependence of the structural, electronic, optical, and chemical properties of 
clusters has been the subject of experimental and theoretical studies in both physics and 
chemistry communities for more than a decade [1;]. Recently, the covalent semiconductor 
clusters have received substantial attention because of their expected importance in micro¬ 
electronics industry 0H3- Since silicon is technologically the most important semiconductor 


material, these studies have attempted to understand the evolution of its structural and op¬ 
tical properties as a function of cluster size. 

Because silicon clusters are highly reactive, they are mostly synthesized in a molecular 
beam under high vacuum conditions |^,[|. Consequently, the structures of these clusters 
are unknown experimentally. Theoretical calculations have now established the structures 
of Sijv clusters in the N = 2-10 atom size range [|l3j. In order to extend this work, we have 
recently carried out tight-binding molecular dynamic (TB-MD) simulations on Sijv clusters 
in the N = 11-14 atom size range. The results of these simulations are reported here. 

The choice of the tight-binding method for the study of the cluster structures is moti¬ 
vated by its accuracy and computational efficiency. Tomanek and Schliiter have shown that a 
properly constructed tight-binding Hamiltonian can yield cluster structures as accurately as 
the first principles methods Jl4|]. More recently, Menon and Subbaswamy have constructed 
a highly accurate tight-binding Hamiltonian for silicon clusters ||H)||TE|| . This Hamiltonian 
includes Harrison’s universal parameters appropriate for the description of bulk Si [|17]], sup¬ 
plemented with two to four additional parameters for the description of the silicon clusters. 
These additional parameters are derived by fitting to the SB bond length and vibrational 
frequency and to the overall size-dependent cohesive energy curve of clusters []T5|,p~6| . Most 


importantly, none of the parameters of this Hamiltonian are fit to any of the cluster struc¬ 
tures. Full details of the Hamiltonian are described elsewhere HT511T6 


As the size of the cluster grows, the number of structural isomers increases exponentially, 
with the result that searching the complete configuration space for the global potential 
energy minimum becomes a formidable task. We found that by combining the tight-binding 
method with the molecular dynamics simulated annealing technique we can efficiently search 
the cluster conhguraton space and determine the ground state geometry. We used this TB- 
MD method in all our calculations reported here. 

Figures 1-4 display the lowest energy structures obtained using this TB-MD technique for 
clusters in the N = 11-14 atom size range. The Sin structure displayed in Fig. 1 consists of 
two tetragons in the anti-prism geometry and three caps. Two of these caps are attached to 
the opposite faces of the two tetragons, while the third one is attached to the edge of the top 
tetragon. Rohlhng and Raghavachari found a similar structure to be a possible candidate 

m . The edge cap in our structure is replaced by a 

0 


We investigated both 


structure for the ground state of Sin 

second face cap in the structure of Rohlhng and Raghavachari 
these structures and found that our structure is more stable by 0.27 eV than the Rohlhng 
and Raghavachari structure. Finally, if we remove the edge cap from the top and relax the 
resulting structure we obtain one of the low energy Sii 0 structures. Hence, we may view 
Sin as a continuation of the Sii 0 structure. Similarly, if we add an edge cap to the bottom 
tetragon we would obtain the twelve atom icosahedral cage structure. Consequently, we may 
also view the Sin structure as an incomplete icosahedral cage structure. 

In Fig. 2 we display the structure of Sin structure. This structure is an icosahedral 
cage. Such a spherical cage structure has not been predicted or observed for any of the 
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12-atom elemental clusters. Adding a face cap to Sii 2 gives the lowest energy Sii 3 structure, 
as shown in Fig. 3. An alternative structure, derived by placing a Si atom inside the cage of 
Sii 2 , is a high energy local minimum. Similarly, adding a cap to the Sii 3 structure does not 
yield the lowest energy Sii 4 structure. As displayed in Fig. 4, Si i4 assumes a layer structure 
consisting of three planes of four atoms each and two adjacent face caps. By suitably rotating 
this structure, we may also describe it as a pentagon sandwich (or prism) with two caps each 
at the top and the bottom. The pentagon prism is somewhat distorted and displaced. We 
also considered a bi-capped hexagonal anti-prism as a candidate for the ground state of Sii 4 . 
However, this structure proved to be unstable, indicating that six-atom ring structures are 
still not favored in these small clusters. 

The Ar^r clusters in this size regime assume different structural pattern than the corre¬ 
sponding SW clusters |I9|. For example, Ari 3 is an icosahedron but with an atom inside 


the cage |T9|- Furthermore, if the central Ar atom is removed the resulting Ar 12 structure is 


unstable. The stable Ar 12 and Ar n cluster structures may be derived from Ar 13 by removing 


one and two surface atoms, respectively (T9 


derived from Ar 33 by adding one and two face caps, respectively, to its surface [19 


Similarly, Ar 14 and Ar 15 structures may be 

Hence, 


we may view the structures of Ar^ (N = 9-15) clusters in the neighborhood of Ar 43 as being 
derived from Ar 13 by the addition or subtraction of surface atoms. Such a simple structural 
relationship exists for Sijv clusters only in a very limited range of cluster sizes. For example, 
Sin and Sii 3 cluster structures may be derived from the Sii 2 structure by the removal or 
addition of one surface atom. However, Siio and Sin structures are not so simply related to 
the structure of Sin- Similarly, surface caps stabilize Spy structures much more than the in¬ 
terior caps. In sharp contrast, removing the interior atom from the Arjy (N = 9-15) clusters 
make the resultant structures unstable. Hence, the Spy structures are quite different from 
those of the corresponding Ar^y clusters. Since, two-body interaction is the dominant part 
of the Ar-Ar potential, the differences between Ar^y and Sijy cluster structures are due to 
the many-body interactions of the Si-Si potential and the directional nature of the covalent 
bonding. The directional bonds in Si^y stabilize the cage structure, while the isotropic Ar-Ar 
pair potentials stabilize closed shell structures with maximum coordination. 

Previous work on silicon clusters in this size regime employed both classical potentials 
|20H22| as well as first principles methods |T3|, |T8|j23j . Our calculated cluster structures in 
the N = 2-10 atom size range agree with the first principles calculations of Raghavachari 
0. Furthermore, our structure for Sin is similar to the one obtained by Rohlhng and 


Raghavachari [|TBJ. However, our structure for Si 43 is qualitatively different from that ob¬ 
tained by Rothlisberger and Andreoni |23|]. The structures of Sii 2 and Si i4 have not been 


investigated using the first principles methods. 

Classical calculations exist on the structures of all these clusters. The calculations of 
Barojas and Levesque ^0[ and Feuston and co-workers [^IJ employed the Stillinger-Weber 
potential [|24|. However, this potential almost always yielded incorrect cluster structures 
because it was derived by fitting to bulk Si properties. Indeed, the predicted structures 
of Sin and Sii 3 using this potential disagree with our results as well as with those of the 
first-principles methods. Chclikowsky also reported the structures of Sii 3 and Sii 4 structures 
using his own classical potential that was designed specifically for clusters |22l]. However, his 


prediction of atom-centered icosahedron for Sii 3 is in disagreement with our calculations as 
well as with that of Rdthlisberger and Andreoni. Likewise, Chclikowsky’s structure for Sii 4 
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is quite different from our structure reported here. There are no first principles calculations 
on Sii 4 - This comparison shows that our tight-binding method gives results consistently 
different from those of classical potentials, while agreeing quite well with the calculations of 
first-principles methods. 

There is a concern that the tight-binding Hamiltonian we employ may be behaving more 
like an isotropic pairwise additive potential. If that were the case, we would have obtained 
atom centered icosahedral structure for Sii 3 , as is the case for Ari 3 flHIfl . However, it turns out 
that the atom centered icosahedron is a high energy local minimum for Sii 3 - Furthermore, 
our calculated structures for clusters in the N = 2-10 atom size range are in good agreement 
with the first principles calculations [|13j. These two observations confirm that our tight- 
binding Hamiltonian correctly incorporates the many-body interactions necessary for the 
accurate description of the cluster structures. 

In summary, our tight-binding molecular dynamic simulations predict unforeseen struc¬ 
tures in the N = 11-14 atom size range. The Sin cluster is an incomplete icosahedron, Si^ 
is an icosahedral cage, Sii 3 is a surface capped icosahedron, and Sii 4 is a bi-capped 4-4-4 
layer structure. Since most of these are low symmetry structures, they may be described 
differently by viewing them through different orientations. All these clusters have low energy 
isomers that may play an important role in the observable properties. While Sin_i 3 are close 
to being spherical, Si i4 is a prolate ellipsoid. The Si i2 cluster in particular assumes highly 
spherical icosahedral cage structure. Such a spherical cage structure has not been predicted 
or observed for any of the silicon clusters. The most characteristic feature of all these clusters 
is that they are all surface. 

This research is supported by the New York University and the Donors of The Petroleum 
Research Fund (ACS-PRF Grant # 26488-G), administered by the American Chemical So¬ 
ciety. 
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FIGURES 


FIG. 1. The lowest energy structure of Sin cluster. It is an incomplete icosahedron. 


FIG. 2. The lowest energy structure of Sii 2 cluster. It is an icosahedron. 


FIG. 3. The lowest energy structure of Sii 3 cluster. It is a surface capped icosahedron. 

FIG. 4. The lowest energy structure of Sii 4 cluster. It is a layer structure consisting of three 
planes of four atoms each and two adjacent face caps. 
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